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Abstract: We systematically study the possibility of determining the spin of new particles 
after their discovery at the LHC. We concentrate on angular correlations in cascade decays. 
Motivated by constraints of electroweak precision tests and the potential of providing a Cold 
Dark Matter candidate, we focus on scenarios of new physics in which some discrete sym- 
metry guarantees the existence of stable neutral particles which escape the detector. More 
specifically, we compare supersymmetry with another generic scenario in which new physics 
particles have the same spin as their Standard Model partners. A survey of possibilities of 
observing spin correlations in a broad range of decay channels is carried out, with interesting 
ones identified. Rather than confining ourselves to one "collider friendly" benchmark point 
(such as SPSla), we describe the parameter region in which any particular decay channel is 
effective. We conduct a more detailed study of chargino's spin determination in the decay 
channel q —* q + — > q + + LSP. A scan over the chargino and neutralino masses is 
performed. We find that as long as the spectrum is not too degenerate the prospects for spin 
determination in this channel are rather good. 



1. Introduction 



Naturalness of the weak scale implies the existence of new physics beyond the Standard 
Model at the scale ~ TeV. Typical new physics scenarios predict the existence of a set of 
new particles at that scale. The Large Hadron Collider (LHC) gives us a great opportunity 
for discovering those particles. 

In order to understand the nature of the new physics, it is necessary to measure its 
properties in detail. One of the obvious tasks is to reconstruct the masses [1, 2, 3, 4, 5, 6, 7, 8] 
and gauge quantum numbers of the new particles from experimental data. A recent study [9] 
demonstrated the challenges of such a goal and suggested possible directions in achieving it. 

On the other hand, there is another, at least equally important, LHC inverse problem: 
how do we determine the spin of any newly discovered particle? The proposed new particles 
in several main candidates of new physics scenarios typically have similar gauge interactions. 
They could often be organized as partners of the known Standard Model particles with the 
same gauge quantum number, such as quark partners, lepton partners and gauge boson 
partners. A typical example is the set of superpartners in supersymmetry, including squarks, 
sleptons, gauginos, and so on. Another interesting scenario is the theory space models inspired 
by [10]. The duplication of the Standard Model states in this scenario comes from introducing 
more copies of the Standard Model gauge group. Typical examples are little Higgs models x . 
Measuring the partners' spin becomes a crucial, sometimes single, way to distinguish those 
scenarios. 

Motivated by electroweak precision constraints and the existence of Cold Dark Mat- 
ter, many new physics scenarios incorporate some discrete symmetry which guarantees the 
existence of a lightest stable neutral particle, LSNP. Well-known examples of such discrete 
symmetries include R-parity in supersymmetry, KK-parity in universal extra-dimension mod- 
els [11], or similarly, T-parity in Little Higgs Models [12, 13, 14, 15]. The existence of such a 
neutral particle at the end of the decay chain results in large missing energy events in which 
new physics particles are produced. This fact helps to separate them from the Standard 
Model background. On the other hand, it also makes the spin measurement more compli- 
cated because it is almost impossible to reconstruct the momentum, and therefore the rest 
frame, of the decaying particles. 

The question of spin determination has been revisited recently. The total cross section 
might serve as an initial hint to the spin of the new particles discovered [16]. This is not 
entirely satisfactory because certain model dependence is inevitable when using the rate in- 
formation. For example, a fermion can be faked by two closely degenerate scalars. Moreover, 
such a determination is only possible if we could measure the masses of the particle using 
kinematical information. As demonstrated in [15, 17], typical "transverse" kinematical ob- 
servables are not sensitive to the absolute mass of particles. One can only deduce the mass 
difference between the decaying particle and the neutral particle escaping the detector. With 

1 Another well-known example is the the extra-dimensional setup with the corresponding KK particles. As 
we learned in the past few years, this is related to the theory space models via deconstruction 
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some assumptions regarding the underlaying model there are more subtle kinematical observ- 
ables which, in combination with the rate information, could determine the spin [17]. To 
what extent this could be generalized to a broader classes of new physics particles is currently 
under investigation. 

Therefore, it is important to investigate other possible ways of directly measuring the spin 
of new particles. The typical way of measuring the spin of a decaying particle starts with 
reconstructing its rest frame from the decay products. Then, the angular distribution in the 
rest frame contains the full spin information, independent of the boost. As discussed above, 
we do not have enough kinematical information to boost to the rest frame of the decaying 
particle if the spectrum contains a LSNP. Therefore, it is natural to consider distributions as 
a function of relativistic invariants constructed out of the decay products of a single decaying 
particle. We will focus on this possibility in this paper. 

Various new physics models always have some detailed differences in their spectra. But, 
such differences are very model dependent. Although in principle they could carry interesting 
information, we will focus on spin determination based on Standard Model partners only. 
What we have in mind are two classes of models with almost identical gauge quantum numbers 
and maximal flexibility in their mass spectra. In other words, they would look very similar 
except for their spin content. 

One obvious scenario is low energy supersymmetry, parameterized by the MSSM with a 
conserved R-parity. 

As a contrasting scenario we consider a framework in which all new particles have the 
same spin as their Standard Model counter parts. The existence of a LSNP, is guaranteed by 
the assumption that all the new physics particles are odd under a certain Z2 parity. Special 
cases of this scenario could be the first KK level of UED or T-parity little Higgs. However, 
what we have in mind is a more generic setup and we will not constrain ourselves to any 
special mass or coupling relations imposed by these two scenarios. To emphasize its generic 
nature, we will call it the Same Spin scenario in this paper. We will use symbols with primes 
to label the new particles in the Same Spin scenario. For example, we will use q' to label the 
quark partner, and so on. We will assume the LSNP in this case is a vector and label it as 
A'. 

Typical new physics scenarios have many complicated decay channels. Many kinematical 
distributions can be constructed from them. One of the main goals of this paper is to present 
a systematic survey of the observability of spin correlations in a wide variety of decay chains 
which are generically present in new physics scenarios. We identify interesting decay channels 
to focus on for spin measurements. It is important to notice, as will be clear from our 
discussion, that the usefulness of any particular channel is restricted to a specific range of 
parameters. We describe the kinematical requirements for each of the channels we analyze. 
There is no obvious golden channel. For different points in the parameter space, we will 
generically have different decay channels available. Therefore, we will have to devise different 
strategies depending on mass spectra of underlying models. 

One of the important tools we develop in this paper is a set of simple rules, which 
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summarizes many well-known results concerning spin correlations. Such simple rules allows 
one to gain insight into the angular correlations in decay, without the necessity of going 
through a lengthy calculation. They can be useful in other, potentially more complicated, 
scenarios than the ones considered hereafter. 

Barr [18] investigated a typical supersymmetry cascade involving a squark decay. He 
found that angular correlations exist between the decay products. Barr's method relies on 
the fact that squarks and anti-squarks are produced unevenly in a proton-proton collider. 
There are several follow-up studies along the same lines [19, 20, 21]. References [20] and [19] 
went further and contrasted supersymmetry with the universal extra-dimensions scenario 
[11], [22]. Reference [20] found that with a mass spectrum given by the SPS point la, the 
SUSY model is distinguishable from the UED case. In their study, they assumed the lepton 
can be perfectly correlated with the correct jet. That might be possible if complete kinematic 
information is available, but in practice seems quite difficult. In general jet combinatorics 
must be taken into account. 

One limitation of these investigations is the need for a light leptonic partner. It must 
be lighter than the second lightest neutral gauge boson partner, such as the second lightest 
neutralino (which must be a wino or bino), in supersymmetry. While true in some special 
benchmark models [23], there is no reason to assume this is a generic feature of supersymmetry 
breaking. In fact, it is more generic to assume otherwise, especially if one is driven by the 
problem of naturalness. 

We present a detailed study of spin correlations in the decay chain q — > q + Cf — ► 
q + + N\ in supersymmetry and its counter part q' — > q + W'^ — > q + + A' in 
the Same Spin scenario. Such a decay chain does not require the leptonic partner to be 
lighter than the gauge boson partner and is certainly more generic in parameter space. We 
will assume that the mass splitting between the two lightest states is greater than m\y,z- 
Therefore, the on-shell decay to W/Z always dominates 2 . Our result shows that it is possible 
to observe spin correlation in this decay chain. As a demonstration of the result of our general 
discussion, we map out the parameter region in which this decay channel is useful. 

As part of the analysis we used HERWIG 6.507 [24] which implements a spin correlation 
algorithm. This algorithm was first used for QCD parton showers [25, 26, 27, 28] and later 
extended by Richardson to supersymmetric and top processes [29]. We supplemented the 
code to include the decays of massive gauge bosons (such as KK partners of the gluon) and 
the details are spelled out in appendix B. To the best of the authors' knowledge, HERWIG 
is the only simulator that implements such an algorithm and is therefore suitable for spin 
determination studies. 

The paper is organized as follows: In Section 2 we try to build some intuition by looking 
at the effects of spin on the angular distributions of simple decays. In Section 3, we present a 
survey of spin correlations in various decay channels. Our detailed study of the decay chain 
with qW ± final states is presented in Section 4. There, we take up the task of constructing 

2 If the mass splitting is less than mw,z, sometimes, the off-shell diagram via a squark and slepton can be 
important. Although it is a special case of the mass spectrum, it is certainly worthwhile exploring it further. 
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an observable signal to distinguish SUSY from the Same Spin scenario in the absence of any 
leptonic partners. Finally, in section 5, we comment on possible future directions and present 
our conclusions. 

2. Simple Spin Correlations 

In this section we review some basic angular distributions from simple decays. These distri- 
butions will serve as building blocks in our understanding of the spin correlations in more 
complicated decay chains which we will consider later. 

2.1 Scalar decay 

A scalar does not pick any special direction in space and so its decay is isotropic. It does not 
mean that the existence of scalars spoils any hope for distinguishing them away from phase- 
space. The production of bosons (via a Z° for example) has a different angular distribution 
about the beam axis than that of fermions. This discrepancy can be employed in determining 
the spin of lepton partners (see for example, [30]). However, in our study we will concentrate 
on a single branch in which case it is not possible to distinguish a scalar from phase-space. 

2.2 Fermion decay 

First, we consider the decay of a fermion ipi into another fermion ip2 and a scalar <fi, via an 
interaction of the form 

Vl<HiPl^\ + yR</np2PRipi (2.1) 

Depending on the model, this coupling could be either chiral, yi / Vr, or non-chiral, yi = yR- 
We will see examples of both cases in our study. 

If the coupling in Eq. 2.1 is chiral, ^2 is produced in a chirality eigenstate. If Y>2 is boosted 
then it is in a helicity eigenstate, i.e., polarized. However, ipi is, in general, not polarized and 
therefore the decay is isotropic, even if the coupling (2.1) is chiral and ^2 is boosted. It is 
easy to see how this comes about. If it is a Left handed coupling, yn = 0, then tp2 is mostly a 
right-handed particle, | [). From the transformation of a spinor under a rotation by an angle 
6 we have that, 

IT) -cos (|) IT>+sm(£) 11) 

I I) __ sin (I) |T>+cos(0 |i) 

The angle 9 is defined with respect to ipi polarization axis. Notice that if ipi is left-handed 
polarized, | f), its decay probability is oc sin 2 (|). On the other hand, if it is right-handed 
polarized, | j), its decay probability oc cos 2 (|) . These decay distributions are shown in Fig.(l) 
as a function of cos (6). Unfortunately, ipi itself is normally not polarized and averaging over 
the two process the decay is indeed isotropic. 
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However, if ipi came from the decay of another particle and that vertex was chiral then 
the situation is different. In that case ip\ is polarized and its subsequent decay is governed by 
a non-trivial angular distribution as shown in Fig. (1). Whether the decay involves a helicity 
flip or not determines the sign of the slope. 



|M| 2 




Figure 1: The decay probability for a fermion into a scalar and another fermion of the same helicity 
(solid-black) or opposite helicity (dashed-red) as a function of cos 9. is defined with respect to the 
axis of polarization of the decaying fermion. 

Next, we consider the decay of a fermion into another fermion and a gauge-boson via an 
interaction of the form 

9l4>2^Pl^iA^ + gi&k'fPRilixAp (2.2) 

As before, we consider the case where ip2 is boosted. If the interaction is chiral -02 is in a 
definite helicity state. The fermionic current that couples to is of the form V>q<t^ ipp. 
If the emitted gauge-boson is longitudinally polarized the distributions are the same as the 
decay into a fermion and a scalar. If it transversely polarized it is precisely opposite (i.e. 
same helicity corresponds to sin 2 0/2 and opposite helicity to cos 2 0/2). 

The most important feature of the fermion's decay is the linear dependence of the decay 
probability on cos 9. It is also clear that chiral vertices must be involved in order to observe 
spin correlations (unless the fermion is a Majorana particle, a possibility we discuss below). 

2.3 Gauge-boson decay 

When a gauge-boson decay (2-body), relativity forces the products to be two bosons or two 
fermions. As is well known, when the products are two fermions the angular distribution is 
given by, 

PtranstcOS 0) = i (l + COS 2 0) fl ng(cOS 0) = ~ (l - COS 2 0) (2.3) 

If a gauge boson decays into two scalars via the interaction 

g<fc d» faA", (2.4) 
the angular distribution has the opposite structure, 

^trans (COS 0) = -(l - COS 2 #) P long (cOS0) = COS 2 (2.5) 
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where the subscript on P denotes the initial 
about the polarization axis. The decay of 
the same angular distribution as Eq. (2.5). 
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Figure 2: The decay probability for a gauge-boson into two fermions (left) and two bosons (right) 
for transverse (solid-black) and longitudinal polarization (dashed-red) as a function of cos 9. 

finite mass effects that come into play when the products are not highly boosted. Those tend 
to wash out any angular dependence of the amplitude. Generically these contributions scale 
as m 2 /E 2 . Therefore, as noted before there has to be an appreciable difference between the 
mass of the decaying particle and its products so that m 2 /E 2 < 1/2. 

The contrast with the previous case is clear as the dependence of the amplitude on cos 8 
is quadratic. It is also important to note that the vertex need not be chiral. 

2.4 Higher spin 

By noting that a rotation by 9 of a state of spin j is given by e ld i ay it is easy to see that the 
amplitude for the decay of a particle with spin j is some polynomial of degree 2j, 

P A (cos 9) = a 2j (cos9) 2j + a 2i -i(cos 9) 2j ~ 1 + . . . + a (2.6) 

The coefficients are such that when we sum over all polarizations A we get, 

^P A (cos#) = l (2.7) 

A 

since an unpolarized particle has no preferred direction. In this paper we concentrate on spin 
0,1/2, and 1 and will not consider higher spin. Nonetheless, this is an important issue to 
address. For example, if the partners of the graviton are indeed detected it would be good to 
know whether it is a supersymmetric spin-3/2 object or a Same-Spin spin-2 resonance. 

3. Angular correlations in cascade decays 

In this section, we present a systematic study of spin correlations in a wide variety of cascade 
decay channels. Aside from the matrix element, the kinematics also play a crucial role in the 
observability of spin effects. We lay out the conditions for observing spin correlations in each 
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of the decay channels we discuss. Whether any of the channels is open or not depends on the 
particular mass spectrum. However, it is not unreasonable to expect several such channels 
to be open in a generic model. This is important because the signal from any one channel 
might not be sufficiently strong. In this case we would have to combine the signal from a few 
channels to obtain a high confidence spin determination. 

We focus on a class of specific kinematical observables. It is constructed from the mo- 
menta of two of the observed final state particles. More complicated decay patterns and 
observables consisting of more than two observable particles could also be interesting. 

A generic feature of this type of observables is that spin information of the intermediate 
particle, which has observable decay products on both sides, in the decay chain always man- 
ifest itself as some polynomial structure in the distribution. Indeed, a particle of spin j, if 
polarized, will result in a polynomial of degree 2j. On the other hand, such a method is not 
useful for determining the spin of any particle at the top or bottom of the decay chain. For the 
same reason, very short decay chains such as q — > q + LSP won't contain much information. 

As discussed in Section 2, a key requirement for the existence of any spin correlations is 
for the intermediate state particle to be polarized 3 . A boost invariant way to know whether a 
particle is polarized or not is to study this question in its rest frame using a direction defined 
by its mother particle and the other decay products. We will see examples of such analysis 
in the decay channels we consider below. From the discussion in the previous section, it is 
clear that there are only a few ways for a particle to be polarized in its rest frame, 

1. For Majorana fermions a spin flip results in a different process with different end prod- 
ucts. We must be able to tell those apart (measuring a leptons vs. anti-lepton) . We 
will see this in detail in the discussion to follow. 

2. Dirac fermions must be produced from a (partially) chiral coupling and decay through 
a (partially) chiral interaction. 

3. In general the spectrum of new particles needs not be left-right symmetric. In this case, 
the interaction of these particles is effectively chiral even if the gauge-coupling is vector 
like. A typical example is the asymmetric QCD production of left and right squarks 
when their masses are very different 4 . 

4. For a gauge boson, it must come from the decay of a boosted particle (in the gauge-boson 
rest frame). 

We will see detailed realizations of all of these requirements in various decay channels we 
study in this section. 

We will organize our discussion in terms of different final states. 

3 Strictly speaking, this requirement only applies for on-shell particles. For off-shell particles, the spin 
correlation could have new interesting properties. We will examine the off-shell decay in Section 3.9 
4 We would like to thank M. Peskin for bringing this fact to our attention. 
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3.1 Weak Decay with qP 1 final state 

We will go through the logic of establishing spin correlation this channel in more detail because 
many other channels can be understood following very similar arguments. 

We first consider the supersymmetric case. Here, the decay to qt^ final states will proceed 
through either Dirac or Majorana fermion intermediate states. We consider the case of a Dirac 
fermion first (i.e., chargino intermediate state) and compare it with the corresponding process 
in the Same Spin scenario where the intermediate state is a massive gauge boson W' . We will 
comment on the case of Majorana fermion briefly. For more details, see [18]. 

In the case of a Dirac fermion there are two possible ways for the particle to be polarized. 
If it is off-shell then one of the helicities dominates over the other simply because m 2 /q 2 / 1, 
where q is the fermion 4- momenta (this possibility is also open to a Majorana fermion). This 
might become important if a decay must proceed through an off-shell particle simply because 
no other channel is available. We will not pursue this possibility further, but we comment on 
it in section (3.9). 

The other possibility for a Dirac fermion to be polarized is when both its mother vertex 
and its daughter vertex are at least partially chiral. As an example, consider the decay of 
a squark into a quark, slepton and anti-lepton through a Chargino, as shown in Fig. (3). 
Using the rules we developed in the previous section it is straightforward to understand what 
angular correlations are expected. 

In the rest frame of the Chargino, the decaying squark and outgoing quark define a 
polarization axis. Since the interaction is chiral the Chargino is polarized. Since the second 
vertex is also polarized, we have a polarized fermion decaying into another fermion (lepton) 
and a scalar (slepton). As we saw before, this decay is governed by a first order polynomial 
of cos 8. However, notice that cos 9 is related to the relativistically invariant quantity, t q i, 



where the last equality only holds in the Chargino's rest frame. We can immediately conclude 
that the relativistically invariant amplitude is at most a linear function of t q i with the sign 
given by the explicit details of the couplings. Of course, this can be easily confirmed by an 
explicit computation of the amplitude, as shown in Fig.3). 

In contrast with the supersymmetric case let us consider the decay chain q' — > q + W'^ — > 
q + + v' , where we have assumed that W couples like a Standard Model W . The relevant 
diagram is shown in Fig. (4). If the spectrum is not too degenerate then in the rest frame 
of the W, both the incoming q' and the outgoing q are boosted and mostly left-handed. 
Therefore, the W' longitudinal polarization dominates over the transverse one. Another way 
of seeing the same thing is to note that in the rest frame of W' the fermionic current can be 
written in terms of the gauge-boson polarizations, 
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Figure 3: The decay of a squark through a chargino involves two chiral vertices. As a result the 
lepton's direction is correlated with that of the outgoing quark. On the right we plot the amplitude 
as a function of t q i for rriq — 1000 GeV, mg = 500 GeV and mc, = 300 GeV. The graph is normalized 
to unit area. 

where we have neglected m q . cl and cr are 0(1) coefficients depending on the precise nature 
of the interaction. Notice the suppression of the transverse polarization with respect to the 
longitudinal one by a factor of m q i/E q i in the amplitude. It is clear that when the fermions' 
mass difference is comparable to the W mass, m q i — m q ~ mw', the resulting polarization is 
negligible, since m q i /E q > ~ 1 in the rest frame of the W . 

Since the W' is longitudinally polarized, its subsequent decay into l + v 1 is governed by a 
1 — cos 2 9. Here, 6 is the angle of the outgoing leptons with respect to the axis of polarization 
defined by the quarks. Therefore, the relativistically invariant amplitude squared must be a 
quadratic function of t q i with a negative coefficient in front of the leading power. Notice that 



Figure 4: When q' decays the intermediate W is longitudinally polarized if the incoming q 1 and 
outgoing q arc both boosted in its rest frame. This in turn will result in angular correlations between 
the directions of the quark and the lepton. On the right we plot the amplitude as a function of t q i for 
m q i = 1000 GeV, mw< = 500 GeV and m v i — 300 GeV. The graph is normalized to unit area. 

a gauge-boson does not require the vertices to be chiral. This is important and potentially 
useful in determining the gluon partner's spin. However, it is also more susceptible to mass 
difference effects (see equation (3.2)). In contrast, the fermionic counterpart remains polarized 
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even when rriq is not very different from mg, as long as the coupling is chiral and the outgoing 
quark is boosted. 

Finally, we briefly consider the decay of the squark into ql^ final states via a Majorana 
fermion intermediate state. The relevant diagrams are shown in Fig. (5). For the propagator 
to flip its spin we must place a mass insertion. However, due to the Majorana nature of the 
Neutralino, this corresponds to a different process with different final states than the one 
without a mass insertion. Therefore, the propagator has a definite helicity for each of the 
processes and there are angular correlations between the quark and the lepton. This fact 
was exploited by A. Barr [18] to determine the spin of the Neutralino and further details 
can be found in the reference. One can study q£ + and qt~ distributions to uncover the 
spin information. However, there is a further complication due the Majorana nature of the 
Neutralino. There is always another diagram starting from anti-squark with opposite sign of 
its charge which contributes to the same process but with the opposite helicity structure as 
shown in Fig. (6). As Barr noted, in a proton-proton collider squarks and anti-squarks are 
produced unevenly and therefore the angular correlations are not washed out completely. 



Figure 5: The two possible modes for a decay through a Neutralino. A spin flip requires a mass 
insertion (right), which results in a different process. 



Figure 6: The two conjugate modes starting from an anti-squark for the decay through a Neutralino. 
3.2 Weak Decay with qq final states 

In principle, the decay into qq final states could contain similar spin correlations, since it just 
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replace the leptons in the second stage of the decay chains discussed in the previous section 
with quarks. However, in general we can not determine the charge of the initial jet. Once we 
are forced to average over the two final states shown in Fig. (5), all angular correlations are 
washed out. 

On the other hand, if the decay products of the second decay are a third generation quark 
and quark partner we could, in principle, recover some charge information. It will then be 
possible to extract some spin correlation from such decay chains. The effectiveness of such 
decay channels require further careful studies taking into account the efficiency of identifying 
charge of the third generation quarks. 

3.3 Weak Decay with qW ± final state 

If the charged gauge boson partner is lighter than the leptonic partner then its decay into a 
and LSNP through a non-Abelian vertex is usually the dominant decay mode. This channel is 
shown in Fig. (7) In the supersymmetric case this coupling is at least partially chiral if tan (3 ^ 
1 and the higgsino is not considerably heavier than the gauginos. If m q — m^ » m q , then the 
chargino is at least partially polarized (with respect to the axis defined by the incoming squark 
and outgoing quark in its rest frame). In this case, since the chargino-neutralino-W coupling 
is also in general chiral, correlations between the quark and the outgoing are present. The 
situation is a little more subtle than that since the contributions from a cascade initiated by 
an up-type partner cancel those initiated by a anti-down-type partner. However, due to the 
initial asymmetry between up quarks and anti-down quarks in the incoming PDFs the signal 
is not washed out. This decay exhibits a linear dependence on the variable t q w = (p q +pw) 2 - 
The corresponding distribution for the Same Spin scenario is very different. In the rest 
frame of the W 1 both the incoming q' and the outgoing q are boosted and are mostly left- 
handed or mostly right handed. Hence the W is longitudinally polarized. As a result, this 
decay exhibits a quadratic dependence on t q w with a positive coefficient. 




Figure 7: The weak cascade decay of a quark partner through the non-Abelian vertex in supersym- 
metry (left) and Same-Spin theories (right). 

We will present a detailed study of this channel in Section 4. At this point, we just 
remark that this is a more generic channel comparing with the channel requiring on-shell 
lepton partner in the decay chain. In fact, the existence this decay channel is based on a 
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very minimal set of assumptions about the spectrum, in which only a heavy quark partner, a 
charged gauge boson partner and a LSNP are present. If the spectrum does not even allow 
for this decay chain, we will not be able to extract any information from weak decays. This 
appears to be the most promising channel. 

3.4 Weak Decay with qZ final state 

There is a similar channel with a neutralino as the intermediate particle and a Z° in the final 
state (due to the higgsino-higgsino-Z coupling). This could be a potentially golden channel 
considering the leptonic decay of the Z°. Unfortunately, there are no angular correlations 
since the N^^PxNjZ^ vertex is not even partially chiral. The Same-Spin counterpart is 
slightly ambiguous. If the intermediate particle is a heavy scalar partner of the higgs, there 
are no correlations. However, if the intermediate particle is some heavy Z' this might be the 
easiest channel to discover. As this is not a very generic case we will not pursue it any further. 

There is an additional complication concerning this process. When the Z° decays into 
quarks, this process is experimentally indistinguishable from the previous one we consider 
involving a W ± . However, in most models it is suppressed by a factor of 10 — 50 with respect 
to the chargino channel owing to the higgsino origin of the coupling. Therefore it does not 
present a serious background to it. 

3.5 Weak Decay with qh final state 

The neutralino could also decay into a Higgs and LSP. This is shown in Fig. (8). In the 
supersymmetric case this process is possible because of mixing with the higgsino. Unfortu- 
nately, the hN\N2 vertex is not chiral and no correlation exists between the quark and higgs 
directions. 

In the Same-Spin scenario this process is realized through the higgs coupling to the heavy 
gauge-bosons g'vZ'^A'^h. In this case, a correlation between the higgs and outgoing quark 
exists and follow the same as those for a massive gauge-boson decay into two bosons (the 
amplitude has a quadratic dependence on the variable t q h with a positive coefficient). 

This channel is quite generic and it is important to investigate it further. In certain cases, 
it might be possible to replace the outgoing quark with an outgoing lepton (for example heavy 
slepton production as discussed below in subsection 3.8). In a sense this is an orthogonal 
channel to that considered by Barr [18] as it relies on the existence of heavy leptonic partners 
rather than light ones. 

3.6 Decay of Gluon partner 

In this section, we discuss the decay of the gluon partner into a quark and the quark partner. 
The quark partner subsequently decays into another quark and missing energy. This is shown 
in Fig. (9). This would certainly be the dominant channel of producing new physics particles 
if gluon partners are present in the spectrum. This diagram might prove to be the dominant 
decay mode into missing energy. Unfortunately, neither SUSY nor its Same-Spin counterpart 
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Figure 8: The weak cascade decay of a quark partner through a heavy neutralino into a higgs and 
LSP, in supersymmetry (left) and Same-Spin theories (right). 




Figure 9: The cascade decay of the gluon partner in supersymmetry (left) and Same-Spin theories 
(right). 

have any spin effects present. The super symmetric diagram certainly does not involve any 
correlations between the two outgoing jets owing to the scalar nature of the intermediate 
squark. In contrast the Same-Spin quark is indeed a fermion, however, its coupling to the 
gluon partner is vector like. Therefore, it is unpolarized and its subsequent decay is isotropic. 
In the present work we will not consider this channel any further, leaving a detailed study to 
a future publication. 

If the spectrum is such that the gluon partner must decay into the LNSP via an off-shell 
quark the situation is quite different. We discuss this issue further in subsection 3.9. 

3.7 Strong Decay of Quark Partner 

Next we consider the strong decay of a quark partner. This scenario is slightly specialized as 
it relies on the existence of a squark heavier than the gluino, but it is still generic enough to 
warrant consideration. The relevant diagram is shown in Fig. (10). If such a quark partner 
indeed exist this will be its dominant decay mode. The supersymmetric case still has no 
angular correlations between the outgoing jets owing to our experimental limitations. As 
discussed above, the Majorana nature of the gluino makes it possible to observe correlations 
without having chiral vertices. There are two diagrams, one with a mass insertion and the 
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Figure 10: The strong cascade decay of a heavy quark partner in supersymmetry (left) and Same-Spin 
theories (right). 

other without. The former involves two outgoing quarks and the latter a quark and an 
antiquark. Unfortunately, all that we observe in the lab are two jets and must average 
over the two contributions. Therefore the decay should have no dependence on the variable 

^<3i<32 = (Pqi +P132) • 

The Same-Spin case, however, does posses angular correlations between the outgoing 
jets and is distinguishable from the supersymmetric one. As discussed above, the Same-Spin 
gluon is longitudinally polarized and we expect the decay to be a quadratic function of the 
variable t qiq2 . The biggest challenge such a measurement faces is the signficant background 
due to Standard Model processes. This may not be an insurmountable impasse. The sub- 
sequent decay of q' together with hard cuts on missing energy might reduce the background 
dramatically. This is an important enough channel with a clear enough signature (if isolated) 
to warrant further study which we hope to address in a future publication. 

Notice that the second stage of the decay chain could involve third generation quarks 
and squarks. In this case, since we might recover some charge information, this decay can 
be useful to determine the spin of the gluino. It is not uncommon that the third generation 
squarks are lighter than the first two generation squarks. In particular, RGE running and a 
large third generation Yukawa coupling usually results in a lighter third generation squarks. 
Therefore, if m,- 3 > M- g or just slightly lighter, and 2 > m,q 3 there could be a significant 
enhancement of branching ratio into third generation quark and quark partners. Even in 
the Same Spin scenario, decaying into third generation quark and quark partner could help 
reduce the combinatorial background. 

As mentioned above, in the case of left-right asymmetric spectrum, the two vertexes are 
effectively chiral and that will effect the angular correlations. 

3.8 Decay from leptonic initial states 

Many of the channels discussed above, involved the weak decay of a quark partner. If lepton 
partners are heavy enough, all such channels can be initiated by a lepton partner decay 
instead. The angular correlations are the same, only we replace an outgoing jet with an 
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outgoing lepton. Such a scenario can be realized through the Drell-Yan production of heavy 
lepton partners. 

Such a cascade has several advantages. First, jet combinatorics is not a problem. Second, 
we gain a lot more information because charge and flavor is now available to us. This is 
extremely helpful. For example, in the weak decay with l,W^ final state, there is only one 
channel to consider and no averaging is needed. 

On the other hand, a lepton partner at the beginning of a cascade is harder to come by. 
It could come from the decay of heavy electroweak gauge boson partner, or a Z' coupled to 
leptons. But, that is more model dependent. We also require there to be several states below 
the lepton partner. This could sometimes require special arrangements. For example, in the 
MSSM, we would require a mass hierarchy such as M2 > > M\ > fi. This avenue looks 
promising in certain regions of parameter space. 

3.9 OfT-shell decays 

So far, we have only considered on-shell decay processes. We saw that if the intermediate 
particle is a Dirac fermion the interactions involved must be at least partially chiral. This 
conclusion is modified if that particle is off-shell. Although on-shell decays usually dominate, 
there are special kinematical regions where we are forced to have off-shell decays. For example: 

1. If the quark partner is heavier than the gluon partner, gluon partner will be forced to 
decay through an off-shell quark partner. 

2. The decay of a gauge boson partner to LNSP will be forced to go through an off-shell 
W/Z and lepton/quark partners if the mass splitting is small. The virtual lepton-partner 
and virtual quark channels could be particularly interesting since it brings in new spin 
information about the lepton partner. It could be important over a large mass range of 
the lepton partner since the decay to off-shell W/Z is usually suppressed by mixing. 

To illustrate this point we consider the first example where the gluon partner decays 
through an off-shell quark partner. The relevant diagrams are shown in Fig. (9). The SUSY 
channel obviously has no correlations since the squark is a scalar. However, in the Same-Spin 
scenario correlations are present. The amplitude for this process is 

x-^ o (m 2 , - q 2 )(q 2 - 2m 2 ,)(q 2 - 2m 2 A ,) „ 

^2\M\ 2 oc 2 ' q 2 2 9 m ±1 t qq + f (q 2 ), (3.3) 

where we neglected the trivial denominator. fo(q 2 ) is some complicated function of q 2 , the 
momentum of the internal quark partner, and the masses. It is irrelevant for this discussion. 
Notice that when the quark partner is off-shell q 2 / m 2 , and the coefficient of t qq = (p q + p q ) 2 
is non-zero. This linear dependence of the cross-section on t qq can in principle be distinguished 
from the SUSY case where there is no dependence on t qq . Further study is needed to explore 
the observability of this effect in different models. 

Similar considerations apply for the other case of a gauge-boson partner decay to LNSP 
via an off-shell slepton. 
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4. Determining spin without leptonic partners 



In this section we will explore the decay of a quark partner into a charged weak partner which 
consequently decays into a and missing energy. Let's begin with the supersymmetric case. 
The squark-quark-chargino vertex is given by [31] (we are ignoring the CKM and super-CKM 
matrices as they are quite irrelevant to the following discussion), 

c q gc+ = -92 (uPniUnC! + U 21 C 2 )d + dP R (V u C^ + V 2l C c 2 )u) (4.1) 

where Uij,Vij are the matrices diagonalizing the Chargino's' mass matrix. We are assum- 
ing that the chargino is dominantly a gaugino and ignore the direct quark-squark-higgsino 
couplings. The more important vertex is the chargino-PU + -neutralino coupling, 

C w-cn = 92 W" Ay {0%P L + O^Pr) Cj (4.2) 

where, 

0§ = --LNuVfr + NaVZ (4.3) 

0^ = ^=N* 3 U ]2 + Nt 2 U n 

and Nij are the mixing matrices for the Neutralino. This interaction is usually at least 
partially chiral (when tan(3 / 1). Therefore we expect the amplitude to have some t q w = 
(Pq +Pw) 2 dependence, with a coefficient given by the difference of the couplings in equation 
(4.3). Indeed, in the narrow width approximation q 2 — > we get, 



1 / m|(m| - ml- 2m^j 



IMP oc^ / (4 " 4) ^ + M?,mi) (4.4) 

where fo(S,rrii) is a polynomial given in the appendix. The 3-body phase-space differential 
volume can be written in terms of q 2 and t q w (see for example [32] ) , 

dPS 3 = 1oa \ 2 dq 2 dt qW (4.5) 

with appropriate kinematic boundaries. In the narrow-width approximation the integration 
over q 2 is trivial and simply removes the denominator in equation (4.4) and replaces q 2 — > m^. 
Therefore, the angular correlations in this decay depend on the difference {a 2 L — a 2 R ). 

In Fig. (11) we plot the ratio a 2 R /a 2 L as a function of tan (3 for a few values of ^-parameter 
with Mi = 100 GeV and M 2 = 500 GeV. This ratio is quite different than unity for most 
choices of the parameters. 

There is one additional complication in the SUSY case. As seen from equation (4.1) there 
are two contributions to any process involving the chargino. One comes from the coupling to 
the up squark, while the other comes from the coupling to the anti-down squark. Therefore, 
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Figure 11: The ratio of the right to left couplings a\/a\ as a function of tan/3 for three different 
values of the /z-parameter. M\ — 100 GeV and M 2 = 500 GeV are fixed. 

as shown in details in the appendix these two contributions differ by the sign of the coefficient 
of t q w- Since we cannot distinguish between a jet coming from a down quark and that coming 
from an anti-up quark, we must average over the two contributions. Fortunately, due to the 
composition of the proton there are more up-like squarks produced than down-like squarks. 
Their ratio in production is shown in Fig. (12) as a function of their mass and gluino mass. 
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Figure 12: The ratio of down squarks to up squarks in production as a function of their mass for 
different gluino masses, rrig = 700, 1200, 2000 GeV (black-solid, red-dashed, green-points). 

Let us contrast this with the corresponding process in Same-Spin theories. The inter- 
mediate particle is a spin-1 Same-Spin W^. As we argued before in equation (3.2) it is 
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Figure 13: Theoretical curves (Left graph) for the q-W correlations in the two models (solid, black - 
SUSY, dashed, red - Same-Spin). Cross-section is plotted against the t q i variable. The right plot shows 
the Monte-Carlo simulation. Both graphs are normalized to unit area. 

dominated by the longitudinal mode. Therefore its subsequent decay is dominated by the 
angular distribution of the second equation in (2.5). Therefore the W + boson is preferentially 
collinear or anti-collinear with the jet. We expect a2{t q w) 2 + a>it q w + a o dependence with 
02 > and a% < 0. The computation is quite involved, but the final expression is indeed, 

l-^l 2 = ( g 2_ 1 M 2)2 (Fo(q 2 ,m) + F 1 {q 2 ,m i )t qW +F 2 (q 2 ,m i )(t qW ) 2 ) (4.6) 

where the Fi's are given in the appendix. It is not hard to show that F2(M 2 ,nii) > and 
Fi(M 2 ,mi) < 0. The shape of the resulting cross-section is plotted in Fig. (13) against the 
corresponding SUSY cross-section. The behavior for small t q w is distinctly different than 
the supersymmetric case. The reason is clear. Small values of t q w corresponds to W + being 
collinear with the jet, which is forbidden in the supersymmetric process, but preferred in the 
Same-Spin case. We have also included in Fig. (13) the results of the Monte-Carlo simulation. 



4.1 Experimental Observable - lepton-jet correlation 

Of course, we have to take into account the fact that W + cannot be observed directly and only 
its decay products can be measured. If it decays to quarks and it is possible to distinguish 
these two jets from the rest by reconstructing the W + (as it is on-shell) the signal might still 
be very strong. We include standard cuts to reduce Standard Model background. We also 
consider the contribution from other new physics processes with identical final states. We 
argue that this signal is still a strong candidate for spin determination 5 . In this subsection, 
however, we concentrate on the other option, namely, the leptonic decay of the W^. The 

A detailed study of Standard Model background with sophisticated jet analysis is beyond the scope of this 
paper. We will begin to address this issue in future publications. 
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main challenge we face is that we cannot reconstruct the W + as the neutrino is unobservable. 
To investigate the resulting signal we used the Monte-Carlo event generator HERWIG [24]. 
We implemented spin-correlations for massive spin-1 particles and the details can be found 
in appendix B. For the Same-Spin production matrix elements we used the ones quoted in 
[20] 6 

We expect much of the difference to be washed out once the W + is allowed to decay 
into leptons. The only observables we are left with are the momenta of the outgoing jet and 
that of the lepton. Therefore, we will plot the cross-section as a function of the invariant 
mass t q i = (p q + pi) 2 . Fig. (14) shows the Monte-Carlo simulation results when the lepton is 
correlated with the jet from it own branch. In practice we have no why to tell which jet came 
from which branch and we amend this below. 

While the behavior is indistinguishable at high t q i, it is certainly very different at low 
t q i. We can attribute this to the fact that in the SUSY model, the W + is very unlikely to be 
collinear with the jet, whereas the opposite is true for the Same-Spin model. 



I Same-Spin 

L - SUSY 

4 - ! 




Figure 14: Monte-Carlo simulation for the quark-lepton correlations in the two models (solid, black 
- SUSY, dashed, red - Same-Spin). Normalized cross-section is plotted against the t q i variable. The 
data sets contain ^ 13, 000 events each. 

This plot suggests that the two models are still distinguishable from one another even 
if the W cannot be fully reconstructed. When taking into account the jet combinatorics by 
pairing the lepton with both jets in the event the results do alter. In Fig. (15) we plot the 
cross-section against t q i for events with two jets and one or two leptons (the two-leptons 

"Except for the matrix element M(qq — > q*q*) which was taken from [33] since the one quoted in [20] 
seemed to give production rates which are too large. Strictly speaking, those are inappropriate for the model 
we consider (as they assume degenerate spectrum and only one mass scale, namely 1/R, the compactification 
radius) . However, nothing in our discussion seems to relay heavily on the precise production cross-section and 
we do not expect any major modification to the conclusions below. To simplify the analysis we switched off 
cluster formation, heavy hadron decays and the underlying soft event. 
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case corresponds to both branches decaying into a lepton). The single lepton diagram still 
exhibits the flattening of the cross-section in the SUSY case for low t q \. This is in contrast 
with the rising cross-section for the Same-Spin model. It is possible that close analysis of low 
t q i can pick up this difference once real data is available. There is no obviously observable 
spin dependence in the second case with two leptons in the final state. 




Figure 15: Monte-Carlo simulation for the jet-lepton correlations in the two models (solid, black - 
SUSY, dashed, red - Same-Spin) with no knowledge of the correct pairing. Normalized cross-section is 
plotted against the t q i variable. The graph on the left is for 2 jets and 1 lepton events. The graph on 
the right corresponds to 2 jets 2 opposite sign leptons events. The data sets contain ~ 17, 000 events 
with 1 lepton and ~ 9, 000 events with 2 leptons for each of the models. 



4.2 Experimental Observable - jet-U/ correlation 

In this section we take up the second possibility, namely that the decays into two jets. 
The advantage is the ability of fully reconstructing the four momenta of the W^'s. There are 
disadvantages as well. First and foremost, very naively the Standard Model background is 
significant. Second, since jets are involved, momentum determination involves some amount 
of smearing. This would affect the reconstruction of as well as the angular correlations. 
Third, it is very hard to distinguish between a and a Z° and so when investigating 
background we must consider both. We do not attempt a full Standard Model background 
analysis, however, there are several reasons why it might be possible to reduce such a back- 
ground. First, hard cuts on missing energy can yield a fairly clean sample of beyond the 
standard model physics. Second, we can easily have additional leptons in the process we 
consider. Simply let the quark partner in the other branch (shown in Fig. (17)) to decay into 
a Z' or W'. It is not hard to imagine other possible channels for producing leptons in the 
final state. The problem of background reduction is beyond the scope of this paper and we 
will assume that some non-negligible set of events, containing new physics, can be isolated 
and analyzed. 
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Figure 16: Considering only 4-jet events this is a histogram of the invariant mass tw,ji f° r each of the 
models: SUSY (solid-black), Same-Spin (dashed-red) . The histogram is normalized to unit area. The 
two data sets contained ~ 9, 000 and 2, 000 events, respectively. The normalized error is approximately 
\ZJV~0.04. 

We consider 4-jet events with a typical event topology shown in Fig. (17) together with a 
possible background from another new physics process. In every event we try to reconstruct 
the from two of the 4 jets and then form the invariant mass t q w with the two remaining 
jets. If more than one pairing reconstructs W it is regarded as failure and the event is 
discarded. Our cuts involve > 200 GeV and rj < 4.0. We make no attempt in trying to 
order the jets by the magnitude of their transverse momentum or some more sophisticated 
ordering (This very naive approach yields a significant difference between the two models as 
shown in Fig. (16)). There can certainly be potential improvements on this measurement by 
using more kinematical information of the jets. The linear behavior vs. the the quadratic 
behavior is still visible on top of the background coming from the wrong pairing. 

There could also be background from different processes in the same model of new physics 
with identical final states. One such channel is shown on the right in Fig. (17) for Same- 
Spin theory and Fig. (18) for the SUSY case. We consider all events with 4-jets. Then we 
construct the invariant mass m| ■ = (p^ + pi) 2 for every possible pair. We used HERWIG 
internal algorithm for jet progenitor formation 7 . We include jet smearing effects [34] using 
the ATLAS specs [5]. The results are shown in Fig. (19). The irreducible background does 
not seem to pose a serious concern as the W peak is clearly visible. This should come as 

7 This object is not a jet, but would become one after a proper experimental jet algorithm is applied. 
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Figure 17: 4-jet event topology for the Same-Spin theory. The diagram we are interested in (left) 
together with a possible Same-Spin irreducible background (right) 




Figure 18: 4-jet event topology for SUSY. The diagram we are interested in (left) together with a 
possible SUSY irreducible background (right) 

no surprise. With a random pairing of such energetic jets, the chance of reconstructing a 
quantity with < 100 GeV is fairly low. This can probably be made even sharper with some 
simple cuts on jet energy. For example, the initial jets from the squark decay tend to be more 
energetic than those coming from the W. 

The above discussion is by no means a proof that the reconstruction of the W is an 
easy task. It merely serves to show that it is not hopeless and to encourage further study 
of this possibility involving a proper detector simulator and a more detailed analysis of the 
background. 

4.3 Scanning Mi and M2 

It is important to map out the regions of parameter space in which the different channels 
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Figure 19: Considering only 4-jet events we plot a histogram of the invariant mass m| fc - ( = (pk +pi) 2 
of all possible pairs of jets. The graph is normalized to unit area. 

are useful. In this section we present the results of a scan covering the subspace spanned by 
(Mi,M2). We also set m\v' = M2 and tua' = M\. While there are other parameters which 
effect the results (such as tan/3, rrig, fj,, etc.) we focus on those two parameters, (M±,M2), 
which are the most determinantal to the observability of angular correlations for the channel 
we consider. 

When performing such a scan, we wish to assign a number to quantify the difference be- 
tween, say, the distribution describing the SUSY scenario and that of the Same-Spin scenario. 
There is no unique choice for such a number. Moreover, such an assignment can be problem- 
atic as it can overlook differences that a more careful analysis would pick up. Therefore, the 
results of this section should be understood with the following proviso in mind: the numbers 
assigned for the different points in parameter space carry only relative importance among 
themselves and have very little absolute meaning. They indicate that in certain regions spin 
determination is easier as compared with other regions. 

There is one more point to keep in mind. We will compare the distributions produced 
for the two models by matching their spectrum. This is incorrect. One should match the 
cross-section first as that is the actual experimental observable. Unfortunately, we had poor 
control over the production rates for the Same-Spin scenario However, this is not entirely 
misleading. As part of our results we will consider the observability of each of the models 
against phase-space. Therefore, if nothing else, we are able to tell when spin effects are present 
at all. 

To quantify the difference between the distributions we use the Kolmogorov-Smirnov test 
for goodness-of-fit. In this test, the cumulative distribution functions (CDF) of both data 

8 The production cross-sections given in [20] for the Same-Spin theory have only one adjustable parameter, 
namely the inverse compactification radius. 
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sets are compared. The D-statistics is simply the maximum vertical difference between the 
CDF's. In other words, 

D = sup\Fi{x) - F 2 {x)\ (4.7) 

where F\(x) and F 2 {x) are the CDF's for the two data sets. The p- value assigned to the D- 
statistics is low when the two data sets come from different underlying distributions. There 
are several advantages to the Kolmogorov-Smirnov test. First, it is a non-parametric test so it 
applies to general distributions. Second, it is independent of the way we choose to histogram 
the data. 

In Fig. (20) we plot the p- value as a function of M 2 for three different values of M±. The 
quark partner mass is m^ q i = 1000 GeV, the gluon partner mass is m^y = 1200 GeV. We 
set tan/? = 10 in the supersymmetric case. For every value of the parameters we produced 
100, 000 events out of which only about 10% passed the different cuts. 

The scan matches our expectations. When M 2 — > M\ the LSP and the W are produced 
at rest and there the correlations with the polarization axis defined by the quark partner are 
diminished. Also, when M 2 — ► m,q q i it becomes harder to distinguish the two data sets. This 
is also expected because in the Same-Spin case the quark partner, q', is not at all boosted in 
the rest frame of the W . Therefore, W is not polarized. 

The results are encouraging. While a more sophisticated analysis (better cuts, better 
fitting, more realistic collider simulation, etc.) is certainly warranted, these initial results 
seem to indicate that spin determination is possible. As pointed out before, one should not 
confine the analysis to special benchmark points, but rather attempt a scan over the parameter 
space. In this way, an inclusive strategy, combining several channels, can be devised to cover 
different exclusive regions of the parameter space. 

5. Conclusions and Future Directions 

We have systematically studied the possibility of measuring the spin of new particles in a 
variety of cascade decays. Generally, the existence of a LNSP renders the reconstruction of 
the momenta of new physics particles impossible. Therefore, we focused on distributions of 
relativistically invariant variables. We identified a set of decay channels which are useful for 
spin determination. 

A general lesson of this study is that even though spin correlations are present in a 
variety of decay channels, the viability of any particular channel is always confined to certain 
kinematical regions. As a result, different models of new physics with different spectra tend to 
give very different useful channels for spin determination. Different strategies will have to be 
employed at the LHC. Therefore, we emphasize that it is important not to confine ourselves to 
any particular benchmark model or any particular decay mode. We should instead explore the 
effective range of all possible decay channels and devise techniques for using them efficiently. 
In this sense, our survey of potentially useful decay channels is the initial step of an important 
task for which many detailed studies still need to be done. 
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Figure 20: A plot of the p- value vs M 2 for different values of Mj.. We set mw = M 2 and tua 1 = M\. 
The p- value is lower when the two data sets are more distinguishable. The other parameters were fixed 
at rriq >q i — 1000 GeV,rrig t g' — 1200 GeV. We set tan/3 = 10 in the supersymmetric case. Both data 
sets contained ~ 10, 000 events. Notice that when the mass splitting between Mi and M 2 is small, 
i.e. the left-hand side of the plot, it is harder to distinguish these two scenarios. Supersymmetry 
and Same-Spin scenario become more distinguishable as the mass splitting increases. However, as Mi 
approaches mq %q i the difference is again diminished. 



As part of this program we studied the decay channel q — > q + C' ± — > q + + LSP. A 
scan over the chargino and neutralino masses was performed. We found that as long as the 
spectrum is not too degenerate the prospects for spin determination are rather good. 

One of the most important challenges is to measure the spin of gluon partners. The 
difficulty in such a study is that the decay products carrying the spin information of the 
gluon partner are usually jets, missing the charge information. In addition, such channels 
usually have larger combinatorial background and Standard Model background. A similarly 
challenging task is a direct measurement of the spin of the quark partner. 

We have considered only kinematical variables constructed out of two objects of the decay 
products. It is in principle interesting to study the possibility of using more complicated 
kinematical variables. 
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Notice that decay channels involving leptons are generally more promising than those 
involving only jets. This result stems from the fact that we have charge and flavor information 
from the lepton, which could help us in separating the channels. We have assumed no such 
information from jets. Therefore, we have limited information from channels decaying into 
quarks, except for the third generation. Moreover, losing charge and flavor information from 
the jet significantly increases the combinatorial background, since typical new physics signals 
for the scenarios considered in this paper almost always contain several hard jets. Therefore, 
any potential information about the charge and flavor of the initial parton of the jet will be 
very helpful in spin determination. 

We would like to remark that reducing combinatorial background is very important in 
extracting more information about the underlaying new physics. This lead the authors of [9] 
to conclude that to what extent jet charge could be measured is an important study for LHC 
experiments. 

We have not studied the possibility of measuring spin in the production. This is certainly 
a very important area to be investigated carefully. Barr [30] has studied the measurement 
of the spin of muon partners using angular distribution from pair production. In principle, 
angular distributions of production of other partners should carry similar information. This 
is a subject currently under investigation. One of the main complications is the existence 
of t-channel productions. Such production channels bring in more partial-waves and tend to 
wash out characteristic angular distributions. 

We would also like to emphasize that although we have only compared SUSY with the 
Same-spin scenario, the general rules we have developed in the paper are easily applicable to 
any generic scenario with different spin content of new physics particles. 

Acknowledgments: We would like to thank the following people for useful discussions: 
Nima Arkani-Hamed, Michael Peskin, Gui-Yu Huang, Jesse Thaler and especially Tao Han 
whose encouragement and suggestions were very helpful. 



A. Matrix elements calculations 



In this appendix we present the explicit matrix elements for the decays considered in the text. 
We begin with the Same-Spin decay channel q' — ► q + W — ► q + W + A'. The matrix element 
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is given by, 



iM = 




m 



w 



m 



w 



X e (g P v(-q - P3)a + 9ua(P4 + O) p + 9ap(P3 ~ Pa)u) ^ (j>4 ) £ P (j>3 ) 



We have borrowed the usual Standard Model coupling for the vertices. The probability 
amplitude is then, 

\ M \ 2 = C 2 (t qW ) 2 + Clt qW + C 



pol 

Where the coefficient functions are, 



C2 



1 rn$y + (m\, - q 2 ) 2 + 2m 2 v {hm\ l - q 2 ) 



(q 2 -m 2 v/ ) 2 ■m 2 v m 2 A , 

1 1 

(q 2 - m 2 w ,) 2 m 2 w ,m 2 Al m 2 w 
x (-(m 2 , - rriq)(m w - m A ,){m A , + lOm^m^, + m A , - q 4 ) 
+ m w ,(—m%r — llm w m 2 A , — Wrr^m^ — m A , + ?>rn\ v q 2 
+ lAmiym^q 2 + 2>m 4 A ,q 2 — ?>rri^ v q 4 — 3m 2 A ,q 4 + q 6 
+ m 2 ,(m w — 10m w m A , ~~ 3?77^ + 4m A ,S — q 4 )) 
+ (7772(-3m^ + m^, - q 4 + m 2 w (-Wm 2 A , + 4g 2 )))) 

The last term, cq is too complicated to present and carries little significance. It is not hard 
to show that C2 is always positive for any real choice of q 2 . 

Next we are interested in the corresponding SUSY decay chain shown in Fig. (21). The 
first vertex is given by, 

-g 2 Vu dP R C c u and - g 2 U? 1 6p l u d* (A-l) 
respectively. The NCW + vertex is given by, 

92W^N l ^{Of 3 P L + 0§P R )Cj (A-2) 
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Figure 21: Two different process contributing to squark decay into C + and a quark. The chargino 
consequently decays into W ± and LSP. 

The matrix element for the decay initiated by a down-type squark is, 

M S = glU^ uM^^Pl + O^P R ) \ +M ^ 2 P L v u {p 2 )eM (A-3) 

q 

ad + bMr, 

= U* n u^p^ f q2 _ M Z v u (p2)eM) (A-4) 

c 

where, 

a^Oft and b = g 2 2 0[ x (A-5) 
The squared amplitude is given by, 

(6 2 m^ - a 2 q 2 )(q 2 -ImL + m%) 

* M f = 2,„^-^)" "V + C-t. (A-6) 

Therefore, in the narrow width limit, the slope depends on the difference b 2 — a 2 . The problem 
is that the second diagram contributing to this process (with an up-squark decay) has the 
opposite sign for this coefficient. Notice that, 



dP R C c u = -{C c ) T P R (Fu = -CC T P R (Fu (A-7) 

= -dp R c T (Fu = -&p R c T aFu 

= &P R d c u 

When contracting this operator with the NCW + vertex we get the opposite spin struc- 
ture, 

Mu = g 2 2 V u u^(p 4 )^(O^P L + 0*P R ) { + % P R v d ( P2 )eM (A-8) 



q 2 -Ml 



hi + aMf, 
Vu UffipA)^ 2 _ 2 v u ( P2 )e^(p 3 ) 
q m c 
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Experimentally, we cannot distinguish between up quarks and down quarks. We must there- 
fore average over the two contributions, 



£ \M\ 2 oc (b 2 - a 2 )(f d (U^) 2 - f u (Vn) 2 )t qW + Const. 



(A-9) 



u,d 



fd and f u is the fraction of events with a down-squark or up-squark, respectively. 
B. HERWIG implementation 

In this appendix we review the implementation of matrix elements into HERWIG. This section 
is relevant to any Monte-Carlo program using the S and F functions of Eijk and Kleiss 
[35] to calculate helicity amplitudes. The usage of these functions results in very efficient 
computations. The price the user has to pay is the complexity of the expressions. These 
functions are then used in the spin correlations algorithm devised by Knowles and Collins 
[25, 28]. We hope to provide a brief but fairly self-consistent presentation below. 

HERWIG is an event generator consisting of roughly 5 phases (for a complete description 
of the program consult [24] ) : 

1. Hard process, where the particles in the main 2 — > 2, 3 event are generated, e.g. q q — > 
g g or e+e" -> qq. 

2. The parton shower phase involving the QCD evolution from the collision energy to the 
infrared cutoff. 

3. Decay of heavy unstable particles before hadronization, such as top quark and SUSY 
partners. 

4. Hadronization stage 

5. Decay of unstable hadrons. 

We are mainly concerned with the third step. In order to decay any unstable particle 
(e.g. gluino, squark etc.), one must provide HERWIG with its different decay channels and 
the corresponding matrix elements. To keep track of spin correlations the matrix element 
must include the external polarizations. For example, let's consider the decay of a heavy 
fermion into a fermion and a scalar (top decay into bottom and higgs). The matrix element 
is given by, 



In order to compute this spinor product, HERWIG requires the user to express it in terms of 
the S and F functions of Eijk and Kleiss [35] defined below. This facilitates the algebra, but 
obscures the expression. Let's briefly recall the construction. 



iM = igu(q, A 2 )w(j>, Ai) 



(B-l) 
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One begins by expressing any polarization in terms of massless spinors (for a pedagogical 
review, see[36]). Two basic 4-momenta, ko,ki are chosen such that, 

k -k = 0, k -h = 0, fci ■ fci = -1 (B-2) 

The basic left and right helicities are then defined via, 

UL(k )u L (k ) = P L ft , u R (k ) = ftxULiko) (B-3) 

The helicity spinor for any other momenta (not collinear with ko) is then given by, 



u x (k) = ^j(i )/V2M (B-4) 

A massive spinor is nothing but the linear combination of two massless spinors of opposite 
helicities. It can be written as, 

(f + m) u^ x (k 2 ) 
n(P ' A) " V2p-k 2 

' {S x (k 1 ,k 2 )ux{k 1 )+mu-x(k 2 )) (B-5) 



where p = k\ + k 2 is decomposed into two massless 4- vectors and the S function is defined as, 

Sx(h, k 2 ) = uxih^xih) = (S_ x (k 2 , fci))* = S x (k 2 , h) (B-6) 

Notice that, \Sx(ki, k 2 )\ 2 = m 2 . 

All expressions can therefore be reduced into products of massless spinors, with the 
possibility of having a 7^ matrix sandwiched in between. It proves useful to define the F 
function as well, 

F{k 1 ,X 1 ,p,k 2 ,X 2 ,M) =u Al (fc 2 ) {f + M) ux 2 {k 2 ) (B-7) 

for some p which is not light-like. Matrix elements for many typical processes were already 
implemented in HERWIG by P. Richardson. As an example, the above matrix element (B-l) 
can be written as, 

M = gu(q, X 2 )u(p, Xi) 

= 9 ^/2p-p 2 ~^q = q^ (^2(91^2)^2(91) +m 2 u_x 2 fe)) (j* + mi) u_ Al (p 2 ) 
1 1 



^J2p ■ p 2 y/2q ■ q 2 

x 9 (S-x 2 (Q2,qi)F{qi,X 2 ,p,p 2 , -Ai,mi) + m 2 F(q 2 , -X 2 ,p,p 2 , -Xi,m{)) 
For a massless spin-1 particle the external polarization e fl (k,X) can be expressed as, 



£ "(fc, A) = uxik^uxik^/x/Ak-h (B- 
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where k\ is any light-like momentum not collinear with k. 

The extra polarization of a massive spin-1 particle adds an extra complication to the 
calculation. Since the only massive gauge bosons in the MSSM are the W ± and Z, it is easier 
to simply insert the entire matrix element (e.g. t — ► b + e + + v e rather than t — ► b + W + and 
W + — ► e + + zv e ). In Same-Spin theories, there are many massive gauge bosons around and it 
proves useful to develop the needed formalism and implement it in HERWIG. 

Looking at the massless polarization (B-8) it is easy to guess the form of a massive one 
in terms of spinors, 

e"(p,Ai,A 2 ) = -4r- «(p,Ai)yt;(p,A 2 ) (B-9) 

where, p 2 = m 2 and Ai ;2 are the usual spinor polarizations. This might seem wrong at first, as 
it seems to imply 4 polarizations rather than the required 3, but as we will see in a moment, 
(+,+) and (— ,— ) both correspond to the scalar polarization. First, let's verify that this 
indeed reproduces the correct polarization sum, 

^2 e /"(P> A i> A 2)e£(p, Ai,A 2 )/(8m 2 ) = ^ Tr(u(p, Xi)u(p, Xi)^v{p, \ 2 )v(p, X 2 )-/ l/ )/(8m 2 ) 

= Tr((p + m)^(p - m) lu )/{Sm 2 ) (B-10) 

PfiPv 



-( 



~9iiv + 



(B-ll) 



It is straight forward to show that (B-9) corresponds to the different polarizations directly. 
First note that, 

e M (p,Ai,A 2 ) 



2\/2m 2p ■ pi 

x {S\ 1 {Pi,P2)u\ 1 (pi) +m u_ M (p 2 )) 7 M (^2(^1,^2)^2(^1) - rn u_ X2 (p 2 )) 
u\(pihi*u\(pi) - ux{p 2 )-i^u x {p 2 ) /2y/2m \ 1 = \ = \ 2 



= < 



u\(p 2 hfMU\(pi) /2y/2m 



Ai = A + A 2 



In the rest frame of the particle we can take p\ = z\p\/2 and p 2 
Pi + Pi = ( m i 0) and the massless spinors are given by, 

/0\ /o\ 



u-(pi) = v / 2m 



1 



W 



u +(Pi) = v / 2m 




1 

W 



—z\p\/2, so that 



(B-12) 



With a few lines of arithmetic one can verify that, 

e„(p, -, -) = ((1, 0, 0, m) - (1, 0, 0, -m)) /(2>/2m) = (0, 0, 0, 1/^2) 

^(p, +, +) = ((1, 0, 0, m) - (1, 0, 0, -m)) /(2y/2m) = (0, 0, 0, 1/V2) 

e M (p,-,+) = 2(0,m,im,0)/(2>/2m) = (0,1/^2,^/^,0) 

e M (p, +, -) = 2(0, m, -im, 0)/(2>/2m) = (0, 1/^2, -i/V2, 0) 
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and so we conclude that e /i (-, +) = e£, e^(+, -) = e£ and e^(-, -) = e^(+, +) = It 
is now straight forward, albeit tedious, to express any matrix element involving an external 
massive gauge boson in terms of S and F functions and implement it into HERWIG. 

We begin with the matrix element for a gauge boson decay into a fermion - anti-fermion 
pair, 

iM = iA x u(k,X 2 )^P x v(q,X 4 ) e M (p,Ai,A 2 ) (B-13) 
The corresponding expressions are, 

V2 



A^((Ai,A 1 ),A 3 ,A 4 ) 



V4p • p 2 2k ■ k 2 2q ■ q 2 
x (A_ Al (F(/c2,-A 3 ,/c,pi,Ai,m 2 ) F(pi, Ai, q, q 2 , -A 4 , -m 3 ) + (pi -^2)) 

+ (Ai - -Ai)) 



v/4p • p 2 2/c • & 2 2c/ • q 2 

<5-Ai(Pl,P2) - ^(pi,^) 



^((Ai,-Ai),A3,A4) = 



mi 

x ( A_ Al .F(A; 2 , -A 3 , k,p 2 , Ai,m 2 ) F(pi,Ai,g,g 2 ,-A 4 ,-m 3 ) 
+ (Pi ^P2, Ai -> -Ai) ) 

The factor of \/2 above is to guarantee proper normalization of the longitudinal mode. With 
a little bit of care it is easy to obtain the other two diagrams /(/) — > /(/) + g-b.. To turn a 
fermion into an anti-fermion or vice-versa simply send m — ► — m. The incoming gauge-boson 
polarization becomes an outgoing one by conjugation, 

e,u e* = ^(P' a 2)7m u (p> a i) = A i A 2 «(p, -Ai)7 M u(p, -A 2 ) (B-14) 

So we simply have to send Ai — > — Ai in the expressions above, and an overall minus sign in 
front of the A4((\i, — Ai), A3, A4) amplitude due to the A1A2 factor. 

The last diagram we consider is the non-Abelian vertex including 3 external massive 
spin-1 polarizations. 

Pi. £fc 

- 

Pi- i W \^ 

p 3 . v, (A). 

= ig(g^(P3-P2r + g up (-Pi-p 3 r + g p ^P2 +piD (B-15) 

x e p (px, Ai, A^ e* (p 2 , A 2 , A' 2 ) e*(p 3 , A3, A' 3 ) (B-16) 
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We give the corresponding expression for the first term only. The other terms can be obtained 
by trivial permutations and using (B-14) for the outgoing polarizations. 



iM D e* fl ( P2 ,X 2 ,X , 2 )e:( P3 ,X 3 ,X' 3 )g^(p 3 -p 2 ye p (p 1 ,X 1 ,X' 1 ) 
-A(X 2 ,X' 2 ,Xs,X 3 )B(X 1 ,X' 1 ) 



(B-17) 



16\/2mi 



m 2 m 3 



where, 



A(A 2 ,A 2 ,A 3 ,A 3 )= < 



X 2 — X' 2 



^| (F(q' 3 , X 3 ,p 3 ,q 2 , X 2 ,m 3 )F(q 2 , X 2 ,p 3 ,q' 3 , A' 3 , -m 3 ) 
-F{q' 3 , h,P3, q' 2 , ^2,m 3 )F(q' 2 , X 2 ,p 3 , q' 3 , A' 3 , -m 3 ) 
+ A 3 -»• -A 3 ) 



1 fS(q2,q' 9 ,-\2)-S*(q 2 ,q' 9 , 



and, 



B(X 1 ,X[) = { 



X 2 - -X' 2 (F(q' 3 , X 3 ,p 3 ,q 2 , X 2 ,m 3 )F(q 2 , X 2 ,p 3 ,q' 3 , X' 3 , -m 3 ) 

+(A 3 -> -A 3 ,g 2 <-> q2)) 



Ai = A' x (F(^i,Ai,p 3 -p 2 ,9i,Ai,0) - F(gi,-Ai,p 3 -p 2 ,gi,-Ai,0)) 
^1 — -a x — f{qi, M,P3 -P2,qi,M,v) J 



Here, the momenta is written in terms of two massless momenta pi = qi + q[, with qf = 
(q'j) 2 = and pf = mf. 
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